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Abstract: We compute analytically the two-loop contribution to the correlation function 
of the Lagrangian with a four-sided light-like (or null) Wilson loop in M = 4 super Yang- 
Mills. As a non-trivial test of our result, we reproduce the three-loop value of the cusp 
anomalous dimension upon integration over the insertion point of the Lagrangian. The 
method we used involved calculating a dual scattering amplitude. Moreover, we give a 
simple representation of the loop integrand of the latter in twistor variables. 
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1. Introduction and main results 

Over the last few years there has been remarkable progress in the computation of observ- 
ables in planar J\f = 4 super Yang-Mills (SYM). These observables include the S-matrix, 
correlation functions of local operators, Wilson loops and combinations of these. An inter- 
esting class of observables is the correlation function of Wilson loops with local operators. 
In particular cases such correlators are fixed by the symmetries of the theory [1], but in 
general they contain useful dynamical information. 

In this paper we consider the simplest correlators not fixed by symmetries: the corre- 
lation function of a polygonal light-like (or null) Wilson loop with four edges and a local 
operator, which we take to be the Lagrangian C of the theory [2,3]. Such a correlation 
function has ultraviolet (UV) divergences characteristic of light-like Wilson loops [4,5]. 
In order to obtain a finite observable, we normalize this correlation function by the same 



- 1 - 



correlator without the Lagrangian insertion. The finiteness of the ratio follows from the 
structure of UV divergences of light-like Wilson loops. Indeed our perturbative results 
agree with this expectation. 

Conformal symmetry implies that the overall scaling dimension of this observable is 
fixed. Moreover, it is a non-trivial function of one kinematic cross-ratio only, 
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is the only cross-ratio which can be formed by the locations of the cusps xi,...,4 (subject to 
the conditions x\ 2 = x 23 = x 34 = x 4i = 0)> an d the insertion point of the local operator 
x$. Note that as a consequence of the symmetry under cyclic permutations of ^i,...,4 , F 
satisfies the symmetry property F(x) = F(l/x). 

F also depends on the rank of the gauge group N and the 't Hooft coupling A = g 2 N '. 
We will consider F in the planar limit, where it has the perturbative expansion 

F W = E(^2-) F {L - l \*)- (1-3) 
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The first non-planar corrections can appear at four loops. The tree- level and one- loop 
contributions to F have been computed in [2,6], with the result 
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In this paper, we compute analytically the two-loop contribution. We obtain 
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where the functions 
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are manifestly symmetric (antisymmetric) under x — > 1/x for n even (odd). 

The correlation functions on the l.h.s. of eq. (1.1) can in principle be evaluated in 
configuration space. However, there is also a dual formulation of the same objects in 
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terms of integrals resembling scattering amplitudes [7,8]. In order to see this, one can 
think about loop corrections to the correlation functions being generated by Lagrangian 
insertions. Formally, i.e. neglecting regulator issues, one then has the integrand of an 
(L + l)-loop four-point scattering amplitude, with L integrations to be carried out. Since 
all divergences cancel in the ratio of eq. (1.1), one can argue that F can be also obtained 
from a dual calculation, where both numerator and denominator in eq. (1.1) are replaced by 
four-point on-shell scattering amplitudes, the numerator having an additional Lagrangian 
insertion. A representation of F in terms of scattering amplitude integrals was given in 
ref. [6]. Note that in this dual representation, UV divergences of the Wilson loop are 
transformed into infrared (IR) divergences of scattering amplitudes. Of course, the latter 
cancel in the final result for F. 

The integrals one obtains are those for the scattering of four massless particles, but 
with the unusual feature of involving an operator insertion at the point x§. If desired, one 
can use conformal symmetry of the Wilson loops (i.e. dual conformal symmetry in the 
scattering amplitude picture) to send this point to infinity. 

We have performed the calculation both using dimensional regularization, as well as 
in a mass regularization setup [9]. The quantify F is expected to be scheme independent, 
see [10], and indeed we verified (numerically) that both calculations gave the same finite 
result. We found that the calculation was simpler in the massive regularization, as expected 
based on previous experience with similar integrals [11,12]. 

Over the last couple of years, twistor techniques have been tremendously successful in 
describing scattering amplitudes of N = 4 SYM at weak coupling, both at the level of the 
integrand, see e.g. [13-16], and for obtaining analytic integrated expressions, [11,12,17,18]. 
A natural question is whether such techniques will also be useful to understand correlation 
functions of local operators. As a step towards this, we give a simpler, twistorial representa- 
tion of the integrand at two loops. This is closely related to similar simplifications observed 
when studying the exponentiation of scattering amplitudes [11]. This representation has 
several advantages, as we shall discuss. At the one-loop level, we see that the result can 
be written in terms of a single finite integral. Being a one-loop integral, the latter is of 
course known. However, perhaps the simplest way of obtaining this result is to derive a 
differential equation [17] for this single- variable function, which can be readily solved. At 
two loops, we find a very compact representation in terms of five integrals, each of which 
has no sub divergences. We show that the remaining overall divergence cancels between the 
different terms. We find it likely that the differential equation technique of ref. [17] and 
related twistor-space methods will allow for a simpler evaluation of these integrals in the 
future. 

Finally, another interesting feature of the above correlator is that by integrating over 
the point where the Lagrangian is inserted, we recover the expectation value of the four 
sided Wilson loop, and in particular, from its most divergent part, the light-like cusp 
anomalous dimension r cusp [4]. Denoting the necessary infrared regulator by A, one obtains 
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a formula (schematically) 



/^r£4~ A ^ rc " !plog2A + 0(logA) ' (1 ' 9) 

Regulator subtleties and the precise form of this identity for two different regularizations 
will be discussed in the body of the paper. This allows to extract the cusp anomalous 
dimension from a finite quantity. Conversely, knowing the cusp anomalous dimension inde- 
pendently, we obtain an integral constraint on the result. We check that the perturbative 
results up to two loops as well as the strong coupling result, give rise to the correct value 
of the cusp anomalous dimension. 

The outline of this paper is as follows. In section two we describe how to obtain the 
analytic result at two loops from the integral representation previously found in [6]. In 
section three we give a twistorial representation for the loop integrals and in section four 
we show that by integrating over the insertion point of the local operator we reproduce the 
correct value of the cusp anomalous dimension. We end up with a summary of our results 
and outlook, while several technical points are relegated to the appendices. 



2. Analytic two-loop calculation from loop integrals 
2.1 Expression in terms of loop integrals 

The expression for F in terms of loop integrals has been written out in ref . [6] . Converting 
to more standard conventions for Minkowski-space loop integrals, we have 
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Here the 8 permutations refer to 4 cyclic permutations of the points x±, X2, x$, X4, and to 
the swap iiOi4,i2<->% 

The one- and two- loop integrals are shown in Fig. 1 and Fig. 2, respectively. In 
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Figure 1: One-loop integrals contributing to with corresponding numerators. 
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Figure 2: Two-loop integrals contributing to F' 2 ' with corresponding numerators. 



dimensional regularization with D = 4 — 2e and e < 0, they are defined by 
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and similarly for the remaining integrals. 

Although the individual integrals are divergent, the final answer for F should be finite, 
as discussed in the introduction. We have performed the calculation both using dimensional 
regularization, as well as in a mass regularization setup [9]. We found that the calculation 
was simpler in the massive regularization, as expected based on previous experience with 
similar integrals. Below we outline the steps that allowed us to find an analytic answer 
in the massive regularization. We verified numerically that the calculation in dimensional 
regularization gives the same finite answer. 



2.2 Description of calculation 

The method we used for the calculation is standard and straightforward, and we only briefly 
mention the main steps. We first wrote down Mellin-Barnes representations for all integrals 
(see e.g. appendix A for an example or [12], where a similar computation was done). We 
did this by introducing Mellin-Barnes representations one loop at a time. This has the 
advantage that it is straightforward to automate and gives relatively compact answers. 
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Having obtained an expression in terms of (multiple) Mellin-Barnes integrals, we pro- 
ceeded to extract the divergences as e — > in dimensional regularization, or m 2 — > in the 
massive regularization, respectively. 

In order to simplify the calculation, we used the fact that the answer is a conformally 
invariant function that depends on xf through the variable x only. Performing scaling 
limits such as e.g. X{ — > oo leaves this function invariant, but it does simplify individual 
integrals. In taking several limits, we were able to simplify the expression to a point where 
one could verify the absence of divergences in F^ 2 \x) analytically. 

We arrived at a representation for the finite part of F^ 2 \x) in terms of a number of 
one-fold Mellin-Barnes integrals, and one two-fold one. The only two-fold Mellin-Barnes 
integral we encountered is 

/(x) = / H^ u r 2 K)r 2 (i + ^r(-^)r(i + ^) 

x T(l + z x - z 2 )Y{-\ - z x + z 2 )(V Q (z 2 ) + ie) , (2.4) 

with —1 < Re(z\) < Re{z2) < 0, and is the polygamma function fy n = di n+1 ^ logT(z). 

One can also reduce this integral to a one-fold one, as we explain presently. First, 
one writes the polygamma functions as a derivative of a T function w.r.t. an auxiliary 
parameter, 

r(i + z)v Q {z) = iim J-r(i + z + 8i)- r(z) . (2.5) 

Then one can see that the z 2 integration can be carried out using the first Barnes lemma. 
When doing this, it is useful to introduce another auxiliary parameter 5 2 in order to separate 
the left and right poles of r(l+zi — z 2 )T{— 1— z\+z 2 ) — > T(l+z\— z 2 +5 2 )T{— 1— z\+z 2 +5 2 ), 
and send 5 2 — > afterwards. (And similarly for the term with T(—z 2 )T(z 2 ) — > T(—z 2 + 
8 2 )^{z 2 + 8 2 ).) In this way, we obtain a one-fold Mellin Barnes representation for F^ 2 \ At 
this stage, one can write the answer in terms of a series expansion that can be resummed. 
We give an explicit example in Appendix A. 

In this way, we arrived at our final result for F^ 2 \ which is given in eq. (1.6). 

3. Twistorial representation 

In the previous section we have presented the integrals that give us the observable under 
consideration up to second order and explained how to evaluate them using Mellin-Barnes 
methods. However, the set of Feynman integrals we used in eq. (2.1) has a drawback. 
Although the final answer is finite, each individual Feynman integral has divergences. Only 
when we sum all of the integrals together the dependence on the regulator drops out. On 
the other hand, we expect that a good choice of master integrals can significantly simplify 
calculations involving Feynman integrals. For example, introducing infrared finite master 
integrals into the set of basis of Feynman integrals, has both conceptional and practical 
advantages, see [11, 16, 17]. For example, one may hope that it will be easier in such 
representations to find powerful differential equations [17]. We give an example of this at 
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the one-loop order. Recall that F is finite in four dimensions. We will therefore perform 
all manipulations in this section in four dimensions. When individual IR divergent parts 
of the answer are evaluated, it should go without saying that these should be regulated in 
a consistent way, in particular when using numerator identities. 

3.1 One loop 

At one loop, one can use numerator identities described in [11] to express the one loop 
contribution to F in terms of a finite pentagon integral. This can be done directly in 
the dual space, but it is more convenient to switch to twistor space. Here we will only 
present the necessary conventions and refer the reader to [13, 14, 16] for a more complete 
discussion. For any point x\ in the dual space we associate a line in twistor space e.g. to %2 
we associate a line denoted by (12) and spanned by two points in twistor space Z\ and Z 2 
, similarly to points X5, xq, x-j we associate lines (AB), (CD), (EF). The points x\, ...,X4 
define a null polygon, and so the corresponding lines in twistor space intersect. Therefore, 
it is natural to associate the point x% for i = 1, ...,4 with a line (£— 1 i), where the number 
i — 1 is defined modulo 4. Using the incidence relations that relate the momentum twistor 
space to dual space we can find that e.g. 

^ 2 _ (ABCD) 
(AB)(CD)' 

where (ABCD) = £ e ijk iZ A Z j B Z^Z l D and (AB) = (ABI^), where /«, is the infinity 
twistor. In a similar way we can rewrite our integrand in terms of twistor brackets. More- 
over, as all integrals we deal with are conformally invariant, the dependence on 1^ drops 
out and we are left only with the 4-brackets. The resulting expression for the integrand 
can be subsequently simplified by the following identity 

(AB13)(CT>24) + (AB24)(CT>13) = (12AB)(UCD) + (34AB)(12CT>) 

- (23AB)(41CT>) - (A1AB)(23CD) - (12U)(ABCD) . (3.2) 

Eq. (3.2) is a special case of an identity described in [11]. Following the described steps, 
we end up with a very simple expression for the one loop integrand 
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The chiral pentagon integral that we obtain in eq. (3.3) is presented in figure 3. For generic 
cf 4 it is given by 

yQ-)(u,v) =log«logu + Li 2 (l-'u) + Li 2 (l-u)-C2 (3.4) 



with u = x\ 5 X2a/ (^25 x i4) ; v = x %> x \zl ( x 35 x ii)- ^ n the limit x\ A — > 0, both u and v diverge, 
with v/u = x fixed. We find that 

= - - lim ^ {l \u,v) = - [log 2 x + vr 2 ] , (3.5) 
2 x\ 4 ^q 4 
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Figure 3: One-loop chiral integral contributing to 



in perfect agreement with the one-loop result (1.5). Of course, this integral is so simple 
that it can be evaluated by many methods, e.g. using Feynman parameters, see e.g. [11,18]. 
A much more elegant way way of computing it is based on differential equations [17]. Here 
we wish to mention that the latter are compatible with the limit xf 4 = 0. In other words, 
the differential equations of [17] directly apply to that case, and one obtains 

xd x xd x F^\x) = ^. (3.6) 

This, together with the boundary condition F^\x — > — 1) = 0, which follows from inspec- 
tion of the twistor numerator, leads to the result of eq. (3.5). It would be interesting to 
compute the integrals appearing at higher loop orders in a similar way. 

3.2 Two loops 

At two loops, we would like to find an identity, similar to eq. (3.2), that gives a chiral 
representation of the integrand, similar to eq. (3.3). In order to find such an identity, 
we proposed an ansatz for the integrand in terms of a set chiral integrals, that could 
potentially represent our integrand at two loops, with arbitrary coefficients. In order to 
fix the coefficients we performed quadruple cuts. Having determined all coefficients, we 
then checked the equality of the two integrands analytically, by expanding out the twistor 
four-brackets. The representation we found is 

F{2) = ll 2 (-4 + 8/6-4J c + 8J d + / e ), (3.7) 

8 perm 

where the 8 permutations refer to 4 cyclic permutations of the points Z\, Z2, Z$, Z±, and 
swaps Z\ <R- Z^Zi -H- Z3. Here, I's are the chiral twistor integrals, presented in figure 4 
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Figure 4: Twistorial representation of two-loop integrals contributing to F^ 2 K 
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and are denned as follows 
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3.3 Checking finiteness 

Let us now discuss the finiteness properties of the new representation presented in eq. (3.7). 
I a is the square of the one loop pentagon integral in figure 3, which is finite, and it is easy to 
see that lb is also finite. Moreover, one can see that none of the integrals I c , 1^, I e has one- 
loop sub divergences. However, they do separately have overall double logarithmic infrared 
singularities. This is very similar to the behavior of the logarithm of the four-particle 
amplitude discussed in [11]. Here, the remaining divergence cancels for the particular 
combination perm (— 4/ c + 8/^ + I e ), making the final answer finite. Twistor space also 
makes it easy to demonstrate the above statements about IR properties, and we devote 
the remainder of this subsection to show this. 

In order to investigate the divergences of the integrals in figure 4 we adopt the parametriza- 
tion used in [20]. The dual conformal integrals discussed in this paper have only infrared 
divergences. These arise when the loop momentum becomes collinear with the momentum 
of an external massless particle. These limits can be conveniently parametrized in twistor 
space. For example, consider the integration variable Xq, which corresponds to Z c and Z^. 
One of the collinear limits is described by Z c tending to Z2, while Z^ tends to a generic 
point on the hyperplane spanned by Z\, Z2, Z3, 

Z C ^Z 2 + 0(e) (3.8) 
Z d -»• aiZi + a 2 Z 2 + a 3 Z 3 + 0(e) . (3.9) 

For the present purpose, we also need to parametrize the 0(e) terms. In doing so, we 
must make sure that the parametrization is generic enough. That is, we do not impose 
any additional constrains, which can make some factors vanish faster than in the generic 
case. On the other hand, we need to restrict the integration variable xq to be real. This 
corresponds to considering an integration bitwistor Y IJ = Zc Z^ of the form 

Y = a\ (12) + a 3 (23) + e (a 4 (12) + a 5 (23) + a 6 (34) + a 7 (41)) + 0(e 2 ) . (3.10) 

Here the bracket (ij) denotes a line in twistor space spanned by Z{ and Zj. Note that Y does 
not contain "non-local" terms like (13) or (24). Moreover, due to the fact that we are deal- 
ing with an integrand that is explicitly symmetric under cyclic permutations of Z\, Z± 
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we only have to consider the collinear limit presented in eq. (3.10). Imposing that limit we 
confirmed that vanishes while I c , Id, I e diverge for generic a's as I c , Id, I e ~ e • How- 
ever, the considered combination of integrals behaves as ^ g perm (— 4/ c + 8/^ + I e ) ~ e°. 
Therefore, it is finite in all collinear limits, even though I c , Id, I e separately suffer from 
logarithmic divergences 1 . 

3.4 Observation on two loop result 

We would like to finish this section on chiral representation of integrands by making an 
interesting observation that relates the two loop result found in section 2 to finite integrals 
recently computed in the literature. Recall that sending x§ to infinity via a conformal 
transformation our integrals become (in general non-planar) integrals for the scattering 
of four massless particles. The planar master integrals for such a process are all known, 
see [22] and references therein. Recently [23] computed two finite two-loop master integrals, 

called 1+4- and J_| with double-box topology using twistor methods. It is interesting to 

note that these two functions together with their transforms x — > 1/x are related to our 
two-loop result in the following way, 



As is well-known, 3(^3 could also be written as a finite two-loop integral, see for example 
section 3.5. of [27]. It would be natural to absorb this into the definition of I++ and 

J_l We suspect that one can find a connection between the integrals for our observable 

from section 2 and the 4-point integrals computed in [23], by sending x§ to infinity via a 
conformal transformation. However, further work needs to be done to clarify this point. 

4. Relation to the light-like cusp anomalous dimension 

The Lagrangian insertion procedure would naively imply 



We wrote "naively" because both sides of eq. (4.1) diverge double logarithmically, due to 
soft-collinear divergences. A valid equation can be written down within a given regular- 
ization. However, F is defined in the limit where the regulator tends to zero. Therefore, 
we can at most hope that re-instating a regulator in eq. (4.1) will allow us to compare the 
leading divergence of its l.h.s. and r.h.s.. 

As we will see, this can be successfully done using a massive or dimensional regulator. 
See also for a discussion of regulator-independent quantities [10], calculated in massive and 
dimensional regularization. 

It is very interesting to note that the combination — 4J C + 87d + I e is finite if and only if we sum over 
the four cyclic permutations. 



F^(x) + 2((2)F^(x) + (F^Xx)) 2 




(3.11) 
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4.1 Massive regularization 

Consider the four-particle scattering amplitude dual to the four-cusp Wilson loop, defined 
with a massive regulator, M4 [9]. Based on the structure of infrared divergences of the 
latter, we expect the following equation to hold, 

I m2 [F]:= f^ —, f[ X) - =A^log(M 4 )+0(logm 2 ). (4.2) 

Here we only need the leading infrared divergences of log(M4), which are given by (see 
e.g. [10,21]) 

log(M 4 ) = -- log 2 m 2 r cusp + 0(logm 2 ) , (4.3) 
with the cusp anomalous dimension 



r 

1 cusp 



E(^) 1 rffi P , (4.4) 

L=l v 7 



2 / \ \ 3 

4-\ 



= 2 (^H#) +11C '(^)' +0(A4) ' (4 ' 5) 

One might worry that eq. (4.2) does not make sense, since we did not keep the dependence 
on m 2 in F. However, one can argue that this additional dependence will not affect the 
leading divergence. We caution the reader that the same argument is slightly more subtle 
in dimensional regularization, as discussed in the following subsection. 

We wish to verify the above relation (4.2) using our result for F^ 2 \ In order to do this, 
it is convenient to compute the auxiliary integral 

f^ ui , X :_ - = 210^^ + 0(108^). (4.6) 

J " nti( i (5+™ 2 ) "p 

We will give the derivation of an analogous formula in dimensional regularization in ap- 
pendix A. This formula is derived for \p\ < 1, but can be extended to other values of p 
by analytic continuation. Due to the x — >■ 1/x symmetry of F we can assume < x < 1 
without loss of generality. Then, we write F(x) as a series in x around and use eq. (4.6) 
to perform the integration. 

In doing so, one sees that only constants and logarithmically enhanced terms in F 
contribute to the cusp anomalous dimension. Note however that we do need to keep terms 
like log k (x)x n , to all orders in n? The one- and two-loop calculations are elementary. 
Technical details of the three-loop calculation are given in the appendix. To three loops, 
we find 



I m2 [F] = log 2 m 2 



+ C(logm 2 ). 

(4.7) 



2 An instructive example of this is the function log(x/(l + a:)) log(l + a;), which could have appeared at 
two loops. Although it vanishes as x — > (and as x — > 00, due to the inversion symmetry), it would give a 
contribution of (2 to the cusp anomalous dimension. 
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Taking into account eqs. (4.3) and (4.4), we see that this is in perfect agreement with eq. 
(4.2). 



4.2 Dimensional regularization 

The same calculation can also be done in dimensional regularization. In order to get the 
correct result however, we have to be careful. The calculation of F has been done for D = 4, 
and we should really do the whole calculation using D = 4 — 2e. Can we still recover the 
leading 1/e 2 term correctly? The answer turns out to be yes, but we need to go just a little 
beyond the e = approximation in the calculation of F. In fact, on dimensional grounds 
p{L—i) must, have dimension (L — l)e. It is important to take this into account. As far as 
the leading pole is concerned, one can see that this effectively amounts to multiplying the 
naive answer at L loops by a factor of 1/L 2 . Hence we expect 



where the r.h.s. again follows from the known structure of divergences of light-like Wilson 
loops [4], where we have multiplied the usual L-loop contribution by L 2 , as explained above. 
Then, with the dimensional regularization version of eq. (4.6), derived in Appendix A, 

^ 5^-4^ + 0(0. (4-9) 

m nil 4 e "p 

we can reproduce the correct answer for the cusp anomalous dimension to three loops, see 
eq. (4.4). 

4.3 Strong coupling 

In ref. [2] the following answer for F(x) was found at strong coupling, 

F(x) = X [2(l-x) + (l + x)logx]^ + ... , A»l, (4.10) 
(1 — x) 6 Att 

We also note the value of the cusp anomalous dimension at strong coupling, see ref. [24,25], 

r CU sp = 70 + ... , a>i. (4.ii) 

Let us now verify the relation between F(x) and r cusp , using eq. (4.2). 

We could employ the Mellin transform of F(x), but we find it easier just to use a series 
expansion near x = 0. Recall that in this approach the non- logarithmically enhanced terms 
in F{x) do not play a role. (They are needed however in order for F(x) to be well-defined 
at x — > 1.) 

There is a subtle point in this calculation, which concerns interchanging the expansion 
of F(x) for small x and the space-time integration in eq. (4.6). If one does this naively, 
one obtains sum of the type ^2 n>0 {— l) n - A slightly more careful treatment, to be given 
presently, shows that this can be interpreted as 1/2 = 1/(1 + 1) = 1 — 1 + 1 — 1 + .... 
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Let us first compute the integral over the insertion for some generating functions, for 
which the convergence of the series is clear, and then differentiate w.r.t. certain parameters. 
In fact we will see that letting x — > ax and using a as such a parameter will be sufficient. 
This will give a result valid for a < 1, which we can extend to a — > 1. Indeed, we find that 
with 

x log x x 2 log x 



we can write 

x 



9i(x) 



-2(1 



92{x) 



(1 + x) logx] 



1 



lim 

a->l 



X 



92(x) 



1 



d d 2 
9i{x) + T-g2{ax) + ^-53 (ax) 



(1-x) 3 
Moreover, we have, for a < 1, 

I m 2[gi(ax)} =0 + C(logm 2 ) , 
I m2 [g 2 (ax)] =21og 2 m 2 log(l + a) + C(logm 2 ) , 
I m 2[g 3 (ax)\ =21og 2 m 2 [-a + log(l + a)} + 0(logm 2 ) , 
so that we arrive at 



I 2 



-2(1 - x) - (1 + x)logx] 



1 



log 2 m 2 + C(logm 2 



(4.12) 

(4.13) 

(4.14) 
(4.15) 
(4.16) 

(4.17) 



(1-x) 3 

Comparing to eq. (4.2) and (4.3), we find perfect agreement with the strong coupling value 
of the cusp anomalous dimension given in eq. (4.11). 

5. Summary and outlook 

In this paper we considered the correlation function of a local operator (the Lagrangian) 
with a four-sided null Wilson loop, in planar N = 4 SYM. This is an interesting quantity 
due to several reasons: it is finite; it interpolates between a scattering amplitude/ Wilson 
loop and a correlation function and it is a non-trivial function, not fixed by symmetries, of 
a single cross-ratio. Hence it is an ideal quantity to try to interpolate from weak to strong 
coupling. 

We computed analytically the two loop contribution to the above observable. The 
result has the expected degree of transcendentality and reproduces the correct value of the 
cusp anomalous dimension. Furthermore, we have given a twistorial representation for the 
result, which possesses several advantages. 

There are several open problems. It would be interesting to understand better the limit 



x 



at loop level. At tree-level, or, what is the same, the level of the loop integrands, 
this is related to a forward limit of a NMHV amplitude. It would be interesting if, despite 
regulator issues, the corresponding limit at loop level was related to results for loop-level 
NMHV amplitudes. The relevant six-point NMHV amplitudes are known analytically to 
the two- loop order [12] . More generally, one would like to understand different OPE limits 
of our mixed correlator and understand which quantities/anomalous dimensions can be 
obtained from the answer presented in this paper. Finally, it would be extremely nice 
to guess a recursion relation/expression for this correlator, for instance by using twistor 
techniques along the lines of [26]. 
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A. Integration of F{x) = x p over the insertion point 

We would like to compute the leading divergent term of the following integral 



I D (p) =st I d —^ — 

J m 2 \\ 



x 



p 



= i \y - Xi\ 2 



2+e 



where D = 4 — 2e and recall that x = {x\ 5 x\^x\ z ) / {x\ 5 x\ 5 x\^} . This integral is finite in 
the range —1 < p < 1. Using Feynman parameters we obtain 

g 6 (sy r. , „^ ^ (aicvnaiai)-" 

nt,| 9 -x,P = r'(i + p)P(i- P )(i) I 4 " JM £«'- 1 1 E ,„,( S - II ).,.' 

where we have introduced s = xf 3 ,t = x\^. After performing the Wick rotation, the 
integral over the insertion point y can be readily done, and we obtain 

i d(p) = -77777 7777^7 —r-s 1+p t l p / daS(y a, - 1) - v ' K ' n/ _ . (A.l) 

uyy> r 2 (l+p)r 2 (l-p) J y ^ \ ai a 3 t + a 2 a 4 s) 4 - D/2 

Now we introduce new variables (see e.g. [27]) a\ = r/iCi,^ = — Ci)i a 3 = 7 ?2C2 and 
otA = 772(1 — C2)- The Jacobian is simply 771772 and the delta- function constraint implies 
r/i + r/2 = 1. The integration over 771, 772 can be easily done and we are left with 

7T-T 2 (- e )r(2 + e ) , !_„ r 1 (CiC 2 r((i-Ci)(i-C2))- p 

lD =^-T(-2e)THl+p)T^l-p) S J d ^ 2 (sClC 2 + t(l-Ci)(l-C2))^ 

We can separate the s, t dependence by using the Mellin-Barnes representation 

w^w = fwsi jU ^ r(A + =)r( ~ 3) (A - 2) 

and performing the integration over £1 and £2- We are left with 

The contour of integration has to be chosen such that all the poles of + z) are to the 
left and the poles of — z) are to the right. We see that the contour is 'trapped' between 
the poles of T 2 (l — p + z) and r 2 (— 1+p — z — e). After analytically continuing the contour 
(and thereby picking up a residue), we can take the limit e — > 0. We obtain 

I D =4-2e(p) = S ^^ + 0(e- 1 ) (A.3) 
np e z 

The limit p — > exactly reproduces the divergence of the massless scalar box function, as 
expected. 
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B. Three-loop cusp anomalous dimension from integration over F^ 2 \ 



Here we give details of the evaluation of the integral over the insertion point in eq. (4.2). 
We found it technically useful to rewrite eq. (1.6) in terms of the more general class of 
harmonic polylogarithms [28], 

FW(x) =-l [24C2#-i,-i(aO - 12Q 2 H- lfl (x) + 24C 2 # ,o(z) - 4tf_ 2 , ,o(*) 

O - 

+ 8#_i,_i, o ,o(a0 - 4ff_i,o,o,o(a0 + 12iJ ,o,o,o(^) - 12C 2 #-2(aO 

+ 8( 3 H^(x) - 4( 3 H (x) + 107C 4 ] • (B.l) 

This has the advantage that it is straightforward to make the logarithmic dependence of 
manifest. We have (e.g. using the algorithm implemented in ref. [29]) 



+ log 3 x 



8 I 2 
2 



( x ) = - ^ ± bg 4 x + 12C2 log 2 x - 4C 3 log x + 107(4 



(B.2) 



+ log 2 x[4#_i i _i(x)] 

+ logx [-8fl_ 3 _i(x) - &H_i _ 2 (o;) - 12C 2 i?-i(x) + 4F_ 3 (x)] 

+ [24C 2 i?-i + 8ff- 3 ,-i(x) + 8#_ 2 ,_ 2 (x) + 8fl_i -3(3) + &CaH-i(x) - 8H. 4 (x)} } 

In order to perform the integration over the insertion point, we proceed as follows. First, 
we can generate any logarithm from powers of x by differentiating formula (4.6) w.r.t. p. 
Second, we use the expansions of the harmonic polylogarithms encountered above in power 
series around x = 0. We have (see [28]) 



H- n {x) 

#-1,0,-1 0*0 



1=1 



E 

8=1 
OO 

E 

i=i 
00 



-1)V 

i 

-1)V 



S 2 (i) - £ 



if ,-i,-i(x)=^^p5 1 (i 
i=i 



-1)V 



j=i 
00 



i=l 
00 

E 



1)V 



(B.3) 
(B.4) 
(B.5) 
(B.6) 



Here 5 p (n) = £? =1 V* P - 

Following these steps, we see that the first line of eq. (B.2) gives a contribution of 

(3) 

56^4 to Tcusp- Moreover, the second, third, and forth lines contribute 18^4,-2^4 and 
—6^4, respectively, while the last line does not contribute. Combining these formulas we 
straightforwardly obtain the result quoted in the main text. 
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C. Evaluating Mellin-Barnes integrals. 

In this section, we explain how to evaluate one of Mellin-Barnes integrals that arise in this 
paper. In section 3, we introduced Mellin-Barnes representation of the two loop answer. 
Employing the conformal invariance of the answer, we were left with a number of one- fold 
and one two-fold Mellin-Barnes integrals. The two-fold integral is 

Hx) = J H^ 1 -T 2 (-z 1 )r 2 (i + ,)rU)r(i + , 2 ) 

x r(l + z x - z 2 )T(-l - Zl + z 2 )(* (z 2 ) + 7*0 > (CI) 

with —1 < Re(zi) < Re(z 2 ) < 0, and ^> n is the polygamma function such that Vl/ n = 
c?i™ +1 ^ logT(z). This integral can be reduced to one-fold integral with use of the first 
Barnes lemma with a little twist described in section 3. Taking particular care to separate 
right and left poles we obtain 

I(x) = -\ r +t °° gL x -i-* T {-l - z x )T\- Zl )Y\l + z x ) 

x (6T(1 + Z!)(2rf E + *o(-l - zi) + tf (l + *i)) 

+ r(2 + zi)(vr 2 - 67I - 12 7 i^o(2 + z x ) - 6^§(2 + z x ) - 6*i(2 + zi))) 

(C.2) 

In the rest of this appendix, we will present how to evaluate this integral, as an example. 
I(x) can be expanded in asymptotic series in the limit of x <C 1, by the well known 
procedure of closing the contour. Due to the factor of x~ l ~ Zl we close the contour on the 
left hand side of the complex plane. The first pole of the integrand from the right gives 
the leading contribution in the limit i<Cl, the second gives the next-to- leading, etc. By 
summing residues corresponding to that series of poles we obtain the value of the integral 
for any value of x. 

In order to extract those residues, we first simplify the integrand in eq. (C.2) with use 
of simple gamma function identities 

= ^ 6(1 + zi) l j (6(2 7 .(l + zl) - 1)* (1 + zx) - »o(-l - *) 

(1 + zl) (fry! - vr 2 + 6tf§(2 + Zl ) + 6*1 (2 + zx)j\. (C.3) 

Integrand in eq. (C.3) has poles at z\ £ Z and involves polygamma functions with the 
following pole structure for z — > n = 0,— 1, — 2, ... 

lim *o(*i) = - , 1 v + (Sx(-n) - te) + (s 2 (-n) + C(2)) (zi - n) 

Zi-tn [zx — n) \ I \ / 

+ (ft(-n) - C(3)) (zx - nf + (s A (-n) + C(4)) ( Zl - nf + 0{(z x - nf) 
lim ^x(zx) = ( 1 v 2 + (S 2 (-n) + C(2)) + 2(s 3 (-n) - C(3)) (*i - n) 

+ s(S4(-n) + C(4)) (zx - nf + 0((zi - nf) . (C.4) 
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With use of eq. (C.4) we found that the residue of the integrand at z\ = — 1 — i for i > 1 is 

Residue(zi = -1 - k) = (C.5) 

k (2_ _ log(x) Ag^+log 3 ^) _ log 4 (x-) _ log 3 toft (fc) 

1 Xj v fc 4 fc 3 3fc 2 6fc 12 3 

2 2 

+ log 2 (x)(s 2 (k) - S M (fc) - -21og(a0(s 3 (fc) - S 2>1 (k) - S 1>2 (k) + ySi(fc)) 
+ 2S 4 (fc) - 25 3)1 (fc) - 2S 2 . 2 (k) - 2S 1 . 3 (fc) - 7r 2 S M (fc) + ^ + 2S 1 (fc)C(3) - — 

In the above expression we use so-called S-series [30] i.e. S p (n) = Y17=l an< ^ Sp,r( n ) = 
SILi S{r)/i p , etc. All nested S-series i.e. S Ptr (n) in the above expression originate from a 
product of two S-series, which we simplify with Sj(n)Sk(n) = Sj^(n) + Skj(n) — Sk+j(n). 
Now, we would like to sum these residues. This requires series expansions of harmonic 
poly logarithms 3 , already briefly mentioned in eq. (B.3), 

= {H. a {x) fora>0 
f-i i a for a = 

i=l k 1+2; 

-H_ a _ b (x) + H- a - b (x) for a > 

iT ^#_f,(x) fora = 

£■-0,-6,-0(2;) - if_ a _b_ c (x) - i3_ _6 ,- c (x) + #-o-6-c(x) for a > 
- r LJT_ 6 ,_ c (x) + T ^-fl'_ 6 _ c (x) fora = 0. 



OO 



(- 


-l) i+ v 




jo 


(- 






jo 


(- 


l) i+ V 


jo 



E L 4r JL ^) 



Summing the expression in eq. (C.5) from i = 1 to i = 00 using eq. (C.6) and adding 
the residue at Z\ = —1, we found that the integral in eq. (C.l) in terms of harmonic 
polylogarithms is 

XlOg 4 (x) (x- l)k)g 3 (x) log 2 (x) / 2 
i (x) = r ; : H-i(x) H ; r Z7T X — 6ii_l _l(x 

w 12 (1 + x) 6(1 + a?) 6(l + x)V ' 



log(x) 
6 (1 + x) 

2(1 - x)#_ 4 (x) 2F_ 3i _i(x) 2if_ 2 ,- 2 (x) 2fl"_i,_ 3 (s) 



6(1 - x)H_ 3 {x) + (x - 3)vr 2 if- 1 (x) - 12fr_ 2| _ 1 (s) - 12fT_i,_ 2 (x)) 



+ 



1 + x 1 + x 1 + x 1 + x 

vr 2 x#_ 2 (x) 7r 2 i2_i_i(x) 2C(3)fl_i(x) ^ 4 (-ll + 17x) 



3(1 + x) 1 + x 1 + x 180(1 + x) 
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